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ABSTRACT 

We study the effects of rotation on the torsional modes of oscillating relativistic stars 
with a solid crust. Earlier works in Newtonian theory provided estimates of the rota- 
tional corrections for the torsional modes and suggested that they should become CFS 
unstable, even for quite low rotation rates. In this work, we study the effect of rotation 
in the context of general relativity using elasticity theory and in the slow-rotation ap- 
proximation. We find that the Newtonian picture does not change considerably. The 
inclusion of relativistic effects leads only to quantitative corrections. The degeneracy 
of modes for different values of m is removed, and modes with £ = m are shifted to- 
wards zero frequencies and become secularly unstable at stellar rotational frequencies 
- 20-30 Hz. 
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stars: rotation 



1 INTRODUCTION 

Neutron stars are objects of extremely rich internal structure. Although their interior structure is still very uncertain, it seems 
that observations and theoretical studies of neutron stars are quite in agreement concerning the structure of their exterior 
parts. More specifically, there is agreement that neutron stars, 1-2 minutes after their formation, are cold enough to solidify 
their exteriors and form a crystal crust thanks to Coulomb forces between the various atomic nuclei. The crust is covered 
by a very thin fluid ocean, while the interior is formed by a super-fluid mantle (up to 5km in size). The composition of 
the core is highly uncertain. The crystal crust extends from the neutron star's atmosphere 1km down where densities reach 
nuclear densities around 1.3 — 2.4 x lO^'' gr/cm'^. The Coulomb forces of the crystal ions forming the crust can be described 
via the shear modulus fi which is inversely proportional to the 4th power of the ion spacing. Since the restoring force is the 
Coulomb force, the periods of the torsional modes will strongly depend on the shear modulus, and its values will characterize 
the spectrum. Up to now , there are only two detailed calculations of the c rust equation of state (Negele fc Vautherin 19731 : 
Douchin fc Haensel 2001 ) . Both were based on an ap proximate formulation ( Strohmaver et ailllQQl ) leading to a typical value 



)pr 

for the shear modulus of ^ = 1 x 10'^" erg cm~'^pi4 ((schumaker fc Thorr 



1 19831). 



In the Newtonian limit and in the absence of strong magnetic fields, [Hansen fc Cioffi l|l980l ) found that the period of the 
fundamental torsional modes (Pq depends mainly on the radius of the star R, the speed of the shear waves Vs and the angular 
index £ via the following relation 



[£{£+l)]'/\s/R 



(1) 



where Vs = {n/p)^^^, and /x, p are the shear modulus and the density, respectively. Torsional modes, which are axial-type 
oscillations, are thought to be more easily excited during a fracturing of the crust since they only involve oscillations of the 
velocity. The velocity field of torsional oscillations is actually divergence-free without a radial component. Torsional modes 
(t-modes) are labeled as etn, where £ is the angular index, while the index n corresponds to the number of radial nodes in the 
eigenfunctions of the overtones for a specific £. 
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Th e shear and torsion al modes are well st udied in Newtonian theory, see e.g. [Hansen fc CiofBl (ll98C \ 



19831 : iLeins 



Carroll et al 
1994 



McDermot et al. 



(1 19861): Istrohmaverl (Il99l|). while there are only a few studies in general relativity dSchumaker fc Thornel 



Messios et al 



2001 



Sotani et al. 2007a : Samuelsson fc Andersson 2007 : Sotani et all 2007bl ). Relativistic 



2 torsional mode period by roughly 30% for a typical 



effects have been found to increase significantly the fundamental 
M = IAMq, R = 10km model. 

The aim of this work is the study of the ef fect of rot ation on torsional modes. U p to now this is studied only in Newtonian 
theory in a single study by Strohmaver 1 199l[ ). see also Lee fc Strohmaver ( 1996h . That study suggested that the frequency 
of a torsional mode in a rotating star is given by 



a — ao 



(2) 



£{£ + 1) 

where (jq is the fundamental frequency of the torsional mode of a non-rotating star and SI is its rotational frequency. Since the 
typical frequency of the torsional modes of non-rotating stars ctq is of the order of 25-40 Hz, it is obvious t hat they will be subjec t 



to the so called Chandrasekhar-Friedman-Schutz (CFS) instability, even for small rotational frequencies ( Yoshida fc Lee|[2001 



This is a quite interesting result for rotating neutron stars since the instability (together with that of the r-modes) might lead 
to further fracturing and/or melting of the crust. 

Recently, there is increased interest in the study of torsional modes because of the belief that Soft Gamma Repeaters 
(SGRs) could be magnetars experiencing star-quakes that are connected (through the intense magnetic field) to gamma 
ray flare activity. Magnetar star-quakes may be driven by the evolving intense magnetic field which accumulates stress and 
eventually leads to crust fracturin g. There are three SGR events detect ed up to now. The first event was detect ed already in 
1979 from the source SGR 0526-66 jMazets et aI.lll979l : lBarat et al.lll983h . the second in 1998 from SGR 1900-1- 14 jHurlev et al ' 



199£ ) . while the third an d most energetic one was observed in December 2004 fro m the source SGR 1806-20 jrerasawa et al 



20051: (Palmer et al.ll2005l ). Analysis of the tail oscillations, for a full discussion see lSotani et al.l (|2007al ) : IWatts fc Strohmaveil 



1 20061 ) ■ revealed the presence of oscillations from a few tenths of Hz up to ab out 2 kHz. In an att empt to fit the observed 
frequencies to the torsional modes of various EoS for the core and the crust, ISotani et al.l (|2007al ) suggested that neutron 
star models should have stiff equations of state and a mass between 1.6 - 2 Mp). Still these results depend critically on the 



interpretation of the order of each mode. For this, see the discussion in (jSamuelsson fc Anderssonll200' 



r 



In the present work, we derive the p e rturbation equations f o r tors ional oscillations of rotating stars in general relativity 
following the approach by Koiima ( 19921 ) ; Stavridis fc Kokkotaj ( 2005 ) in the Cowling approximation which neglects pertur- 
bations of the space-time. The article is organized as follows. In the next section, we describe the general-relativistic equations 
that have been used to describe the background stellar configuration, and we derive the perturbation equations for torsional 
oscillations of rotating stars in the Cowling approximation. In the third section, we describe the numerical techniques used to 
calculate the modes together with a toy problem. The last is meant to explain the numerical results. The article closes with 
a summary and discussion. 



2 PERTURBATION EQUATIONS FOR TORSIONAL OSCILLATIONS 

We consider a slowly rotating relativistic star described by the metric 



''^"dt^ + e^^dr^ + r^idd'^ + sin^ ddfj)"^) - 2Lur^ sin^ edtd(j) (3) 



-e 



where v, X and uj (the frequency of the local frame) are functions of the radial coordinate r. Up to first order in fi, the 
background 4- velocity of the star is given by 

u"" = [e-'',0,0,ne-"] (4) 

where SI is the angular velocity of the star. The background stellar models are solutions of the TOV equations and an equation 
describing the dragging of inertial frames 

zu" - (^i^' + A' - - 167re^^ {p + e)zj^O (5) 

where zu := ^1 ~ u; . We assume that the star consists of a perfect fluid described by the energy-momentum tensor 

T^v = (p + >^)uf,u„ +pgfii' ■ (6) 

We also assume that the star is isotropic, therefore, the background shear tensor vanishes. 

Due to the spherical symmetry of the background, the perturbations of the backgr ound configurat i on can be decomposed 
into s pherical harmonics. This leads to a large system of partial differential equations Koiimal ( 1992[ ): Stavridis fc Kokkotaj 
( 20051 ). Here the space-time perturbations are omitted (Cowling approximation) and also the coupling between spheroidal 
(polar) and toroidal (axial) perturbations since they only marginally affect the eigenfrequencies of the torsional modes. 
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Under these approximations, the radial component of the perturbed velocity field and the variations of pressure and 
density will remain unaffected since they are polar perturbations. The perturbation of the 4-velocity Su" is related to the 
displacement vector i^" through the relation 5u°' = CuC,"^ ■ For an observer co-rotating with the star, this is translated to: 



5u'>' 
5u 



sin d(j> 
1 dY,^ 



(n — Lj)r e ' Z sin 9 



mrr. 

de 



(7) 

(8) 
(9) 



In other words, the velocity perturbations are described by the time derivative a function Z = Z{t, r) related to the displace- 
ment vector here the dot stands for the temporal derivative. The perturbed energy- momentum tensor, which includes the 
contribution from shear, is given by 



ST^u = {p + (-){5u^Uu -f u^Suu) — 2fi5S^v 



(10) 



where is the shear tens or, Su^ is defined by (Jm^ = CuS^v. Here cr^^ is the rate of shear given by the Lie derivative of the 
shear along the world lines (jCarter &: Quintana 1972 ') 



p°' -L „ p° 



1 



and P^tu is the projection tensor 



(11) 



(12) 



The speed of shear waves on the crust depends on the shear modulus /i, the density e and the pressure p of the star according 
to 111 = ^i/(p + e). A typical value of the speed of shear waves is Vs ~ 10*cm/s. The perturbation equation for the energy- 



momentum tensor 5{T^^ 



provides a single equation for the axial perturbations 



'7" I I ^ I ' \' I 



r^l 2 

z — e 



,A-2, 



1 M 2 



('-!) 



where A = +1)- In the absence of shear = 0, «s = O) , an inertial observer obtains 

X 



-m I Q 



(13) 



(14) 



This gives the well known relati on for the r-mode frequency in the Ne wtonian limit (fi = zu), while in the relativistic case it 
leads to a continuous spectrum (jKoiima Il998l : iBever fc Kokkotaslll999l ^ . In the last case couplings between higher I and polar 
perturbations need to be included for a proper study of the spectrum. 



3 ESTIMATES OF THE MODE FREQUENCIES 

In this section, we present an approximate solution of the boundary value problem and describe the numerical techniques 
used in the calculation of the frequencies of torsional modes of rotating stars. 

For the numerical estimation of the frequencies, we use two different techniques. The first approach assumes a harmonic 
time-dependence of the perturbations. This leads on a boundary value problem. The second approach is based on a direct 
time evolution of equation (|13p followed by a Fourier transform in time of the obtained values at a fixed radial position. We 
will describe only the first approach. The second approach has only been used for the verification of the results. 

The Fourier transform of equation (jlSp ,i.e., the assumption that Z — iaZ, leads to the following differential equation 

,7// , / 4 , , u'\ , 2\-2v J C^r r, CTr r 1 2 /, 2 \1 2,y A - 2 1 

^ +(,7 + '^ +7rj^ |^-2'-i;:^[A+"4'-AjJ^^ x^|^ = °- 

The boundary conditions in the center (or at the lower end of the crust) and on the stellar surface are 

Z ~ /"^ for r -> or = for r Rc and Z' ^ for r^R (16) 

where Rc is the distance of the lower end of the crust from the center. 

The above system of equations defines an eigenvalue problem for the frequencies of the torsional modes a. For its solution, 
we use two approaches. The first is an approximate analytic solution, and the second is a numerical solution. 
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3.1 Approximate solution 



Here we provide an approximate analytic solution to the eigenvalue problem based on Bessel functions. In order to be able 
to use Bessel functions and to treat the problem semi-analytically, we need to simplify equation (I15|) . For this we make the 
following simplifying assumptions, e'^ — ~ 1, ^' — v' — \' ~ Q, vj = Q. and v1 <^ 1. This reduces equation (|15p to a Bessel 
equation 



Z" + -Z' + 
r 

where 







2m Q, Or 

A Vs Vs 



(17) 



(18) 



Together with the boundary conditions (|16p . this equation leads on a simpler boundary value problem. Actually, for high 
frequencies, the eigenvalues could easily be estimated using the WKB method. On the other hand, since the frequency of the 
fundamental torsional mode is relatively small, it is unclear whether the WKB approximation is applicable, in particular in 
the discussion of CFS instability. The last involves the investigation of the limit E ^ 0. 
The general solution of equation (|17|l can be given in the form of Bessel functions 



Z{r) = cir-''^''Ji^i {T.r) + C2r 



-3/2 



(19) 



, the contribution 



where ci and C2 are arbitrary constants. Since the condition for regularity in the center demands that Z ~ 
of the Bessel functions Yn{r) is excluded since it is divergent for r ^ 0. The second boundary condition demands that Z'{R) = 
which leads to the following transcendental equation 



(^-l)J,+ l (fc)-fcJ,+ 3 (fc)=0. 



(20) 



where k — T.- R. The roots of this equation can be found by elementary numerical methods. This leads to the following values 
for k 



1 = 2 

e = 3 



2k„ = 2.501, 7.136, 10.515, 13.772, 
sfcu = 3.865, 8.445, 11.882, 15.175, 
4fcu = 5.095, 9.713, 13.211, 16.544, 



(21) 
(22) 
(23) 



The eigenfrequencies can be obtained from the roots of the equation 



2mD, 



fekn V 



(24) 



by setting fCTn'"' = ik„Vs/R (the frequency in the absence of rotation). In this way, we get an approximate form of the 
torsional mode frequency ar in the rotating frame 



C^n = to., 



(0). 



1 + 



+ 



_fa„(0) ^(^ + l)_ 
which for an inertial observer [a = dr — mQ) has the form 

mfi (i'^ l) 




1 mSl 



(25) 



(26) 



£(£+1) 

Both relations agree with the Newtonian results of Strohmaver (1991), equation (|25p and of Yoshida fc Le3 (2001) equation 
()26|) . Note that eko ~ \/ £{£ + 1) leads to a well known form of the frequency of the fundamental torsional mode for non-rotating 
Newtonian stars, see equation ([TJ. 



3.2 Numerical solution 



The numerical solution of the eigenvalue problem defined by equations (|15p - (|16p is solved by a shooting method. Its results 
ha ve been tested by direc t numerical evolution of the time dependent equation (|13p . The numerical results verify the suggestion 
by Yoshida fc Leel (2001) that the torsional modes are CFS unstable. 

In Figure 1, we plot the frequency of the fundamental torsional mode (2^0 for m = —2. ..2) as a function of the stellar 
rotation frequency Q. We also show the results derived by the approximate analytic method described earlier. It is clear that 
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Q. (Hz) 

Figure 1. The frequency of the torsional mode 2*0 as function of tlie rotation (both in Hz). The dashed line corresponds to the 
frequency of the toy model and the continuous line to the frequency of the relat ivistic model achieved by a direct numerical solution of 
the eigenvalue pro blem. The data are f o r a st ellar model with a realistic EoS A l|Pandharipandelll97ll) for the fluid core and an EoS for 
the crust given by iDouchin fc Haensell l|200ll ') (model A + DH14,). This stellar model has radius R = 9.49 km M = 1.4Mq and crust 
thickness Ar/R 5.15%. 



there is a difFerence in the results of the order of 30% . Mainly, this is due to the fact that in one case we used the approximate 
Newtonian form of the equation, while in the other case we used its exact relativistic form. 

It is obvious th at the torsional m odes become secularly (CFS) unstable even for very slowly rotating relativistic stars. In 
a Newtonian study (jStrohmave il ll99lh torsional mode frequencies, measured in a co-rotating reference frame, can be described 
approximately by the relation where a — ar- Moreover, an inertial observer will measure 



mQ 

e{e + i) 



(27) 



The CFS instability sets in when the frequency of the torsional mode for the inertial observer di = or, equivalently, when 
the phase velocity of the mode is equal to the rotational frequency ar/m = fl, i.e. , when the critical rotation frequency of 
the star is given by 



Q 



inst 



^(^+1) 



(28) 



The approximate results derived earlier (|25|l - p6| l verify the Newtonian results. Also the numerical results agree extremely 
well on the effect of rotation on the frequencies of the torsional modes, i.e.. 



\dn. 



Numerical 



V 0^ J Approximate 



(29) 



with a typical error of the order of 2-5%, depending on the compactness of the star. 

The perturbation equatio ns have been s olved f or a number of different equations of state for the fluid core and the crust 
which are listed in Table 1 of ISotani et al.l (|2007al ). The overall picture is the same, i.e., shows a 20-30% difference of the 
frequencies of the relativistic and the Newtonian equations, and a 2-5% difference in the rate of th e frequency as functio n of 



rotation (^f). The frequencies plotted in Figure 1 are for a stellar model with a re alistic EoS A ( Pandharipandel 1971 1 for 



the fluid core together with an EoS for the crust given bv lDouchin fc Haensell (|200ll ) (model A + DH14) 
has radius R = 9.49 km M = I.4M0 and crust thickness Ar/R f» 5.15%. 



This stellar model 
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4 RESULTS AND DISCUSSION 

In this article, we showed by means of numerical and s emi-analytic me t hods that the torsional modes of rotating relativistic 



stars are subject to the CFS instability as suggested bv lYoshida fc Led (|200ll 'l for Newtonian stars. This instability might be 



only of academic interest since viscosity works against it, and therefore it will probably never prevail. 

We would like to emphasize that the CFS instability does not operate in the up to now observed SGRs because they are 
very slowly rotating stars with periods of the order of seconds. Still it is possible that the torsional modes of the newly born 
neutron stars (soon after they form a crust) will be CFS unstable because rotation periods of the order of tenths or hundreds 
of Hz are expected. The instability of torsional modes in the crust and the unstable r-modes in the fluid core might work 
together in the direction of breaking or melting the crust. If strong magnetic fields are present, then the accumulated stress 
might enhance the above scenario suggesting that the young neutron stars with strong magnetic fields will probably have a 
more frequent flare activity. CFS type rotational instabilities of magnetic field modes will also be present in rotating neutron 
stars, and their effect in the stellar flare activity is part of an extension of this work. 

Finally, rotation (as well as the presence of strong magnetic flelds) will produce frequency shifts towards both lower and 
higher frequencies which poses extra difficulties for the identification of the various observed frequencies from SGRs. 
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